Mediation of Long-Range Attraction Selectively between Negatively- Charged Colloids 

on Surfaces by Solvation 
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We propose a mean-field analytical model to account for the observed asymmetry in the ability to 
form long-range attraction by the negatively charged colloidal particles and not their equivalently 
charged positive counterpart. We conjecture that this asymmetry is due to solvation effects, and we 
phenomenologically capture its physics by considering the relative strength of this water-induced 
short-range repulsion between the different charge species. We then apply our model to the colloidal 
system of negatively charged disks that are neutralized by a sea of counterions and strongly absorbed 
to an interface in a compressible binary system. We demonstrate the resulting coexistence between 
a dilute isotropic ionic phase and a condensed hexagonal lattice phase as a function of density and 
interaction strength. 

PACS numbers: 61.50.Ah, 62.20.Dc 



Colloidal particles have always been the prime con- 
stituents of prototypical self-assembling systems. Their 
phase properties can be easily probed by optical tech- 
niques and readily manipulated via chemical means in 
the laboratory. Their rich chemistry leads to such di- 
verse applications as uses in coating materials, ceramic 
precursors, and the designing and manufacturing of bio- 
logical macromolecules. It has been found that, biologi- 
cally, absorption of colloids into liposomes can be used to 
stabilize vesicles [![. The collective behavior of colloidal 
particles has also been recently demonstrated to provide 
a coorperative probe of membrane surfaces @ • 

Specific colloidal interactions that influence the sys- 
tem's overall properties include Van der Waals disper- 
sion forces, which are the direct consequence of the in- 
duced dipole-dipole interactions created from the quan- 
tum fluctuation in the charge densities electrostatic 
interactions between the charge densities themselves [J, 
as well as an effective attraction of the larger particles 
due to depletion when the colloidal suspension contains 
particles of different sizes @. At the mean-field level 
where we neglect fluctuation effects, it is well known that 
the Poisson-Boltzmann framework precludes long-range 
attraction between like-charges, in contradiction to ex- 
perimental findings @] such as those found in confined 
systems Q • So one may account for these long-range at- 
traction between like charges by going beyond mean-field 
approaches and take in account fluctuation effects. 

However, there has been recent experimental indica- 
tions that this long-range attraction occurs only between 
colloidal particles of negative charge and not of the op- 
posite kind Q. Depletion cannot account for this asym- 
metry since it is entropic in nature and should not dis- 
tinguish between the two kinds of charge. On the other 
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hand, Van der Waals forces are much smaller in magni- 
tude by comparison to be a possible significant contribut- 
ing factor. 

In this paper, we propose a simple, mean-field an- 
alytical model that would demonstrate this asymmet- 
ric behavior between the negatively (positively) charged 
colloidal disks in their ability (inability) to form long- 
range attraction based on their interaction with water 
molecules. In particular, we conjecture that due to the 
bent-core shape and the charge distribution of the water 
molecule, the greatest aggregation of them - swelling - 
occurs between the larger-sized colloidal particles that 
are negatively charged and the smaller-sized counteri- 
ons that are positively charged. Swelling also occurs, 
albeit to a lesser degree, in the opposite case where 
the system consists of equivalently charged positive par- 
ticles with negative counterions. In contrast, swelling 
would occur much less favorably between particles that 
are of the same charge type no matter their respective 
size. Therefore, this interaction with water molecules 
induces an effective repulsion between opposite charges 
and correspondingly a comparative attraction between 
like charges, which exactly counter-balances the effects of 
electrostatics. Consequently, there exists a general com- 
petition between these two types of interaction in the 
system. For a range of interaction parameters and den- 
sities, we will demonstrate the resulting coexistence of a 
low-density isotropic phase and a high-density hexagonal 
lattice phase in our mean-field model. 

As mentioned, the interaction between the colloidal 
particles and water molecules is central to our conjecture 
in accounting for the observed attraction between neg- 
atively charged colloidal particles and the lack thereof 
between positively charged ones. To understand this 
asymmetry, we consider the four scenarios represented 
in Fig. [TJ In the first two situations respectively de- 
picted in Figs. [T]A and QJ3 when we have like-charge 
particles in solution, there exists great frustration in 
the arrangement of water molecules around neighbor- 



ing particles irrespective of the sign on the charge of 
these particles. While it is much more preferable for wa- 
ter molecules to situate themselves between oppositely 
charged particles, there does exist a difference in the 
two cases in which the larger-size particles are negatively 
charged particles (Fig. \Vp) and in which they are posi- 
tively charged (Fig.flp). Due to the bent-core geometry 
of the water molecules and the charge distribution on the 
molecules (positive charge-concentration on the hydrogen 
atoms and negative charge-concentration on the oxygen 
atom), it would be easier for water molecules to orient 
and anchor themselves more energetically favorably be- 
tween a larger particle that is negatively charged and the 
corresponding counterion(Fig. \Vp) than in the reverse 
scenario when the larger particle is positively charged 
(Fig. HP). Thus, we would expect the strongest swelling 
in the scenario depicted in Fig. [TJC and a correspond- 
ingly weaker case in the scenario depicted in Fig. [Tp. If 
we parametrize this water-induced short-range repulsion 
between the charged particles by Wij , for colloidal parti- 
cles of type i = {±} and counterions of type j = {±}, we 
would expect our conjecture to predict the general rela- 
tion {u> ++ ,w__} < < w |_. The detailed elucida- 
tion of their magnitude would necessitate finer knowledge 
in such microscopic details regarding the interaction be- 
tween the colloidal particles and water molecules as their 
precise charge distributions, their associated (quantum) 
fluctuations, and the molecular geometry of the bent- 
core water molecules. Our conjecture is in line with the 
theoretical calculation presented in [8| for the effective 
pair potentials (EPPs) between ions immersed in water 
based on the reference interaction site model (PJSM) 
approach There the EPPs were found after con- 
tracting out the degrees of freedom associated with the 
oxygen and hydrogen sites within the water molecules. 
In their system of NaCl, Gonzalez-Mozuelos found that 

w ++ < w < w |_. Though in his case, the choride 

ions are large than the sodium ions by default, and the 
quantity W-\ was not applicable. 

In what follows, we will present a phenomenological 
minimal model that would capture the selective forma- 
tion of long-range attraction by the negatively charged 
particles in aqueous solution. We will demonstrate that 
our model does indeed yield coexistence between a di- 
lute phase and a condensed phase when the ratio of the 
water- induced repulsion strengths, w±^/w±±, between 
opposite and like charges becomes high enough. For con- 
creteness, we will consider a compressible binary system 
consisting of charged colloidal disks of radius R sub- 
merged in aqueous solution of counterions of radius a. 
Based on our conjecture, when the colloidal disks are 
negatively charged, we would expect the magnitude of 

w |_ > a + ) to be larger than (i?+ > <x_), which 

corresponds to the opposite case in which the colloidal 
disks are positively charged. Again, our model cannot 
predict their magnitudes without incorporating further 
the molecular details regarding the particles' interaction 
with water. Without loss of generality, we will now focus 
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FIG. 1: Interactions between the bent-core water molecules 
(blue) and the positive and negative particles. There ex- 
ists great frustration in the arrangements of water molecules 
around neighboring like-charge particles in (A) and (B). Due 
to the bent-core geometry of the water molecules, they ar- 
range themselves in more favorable configurations of the 
larger negative particles with smaller positively charged par- 
ticles (C) than of the reverse scenario (D). 



on the case when the larger colloidal particles are neg- 
atively charged and demonstrate their formation of an 
ordered phase in our model. 

For the dilute isotropic phase, its homogeneity and 
low density of charged particles permit its treatment as 
a two-dimensional charged gas in a homogeneous back- 
ground that displays nonselective interactions with the 
charged components. As previously mentioned, we have 
two competing interactions in the system: the water- 
induced short-range repulsion between the charged disks 
and the electrostatic interaction between them. Based 
on the EPPs numerically computed in the RISM frame- 
work [1, H|, we analytically model the short-range re- 
pulsion between the counterions and colloidal disks as a 
Gaussian curve enveloping the numerically plotted func- 
tional forms of their EPPs [8|, |9( , which in Fourier space 
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would take the form 

U s k r = 47TR 2 W—e- k2R2 +4TTa 2 w ++ e~ k2R2 

+8TraRw-+e' k2aR , (1) 

in units of fcgT, The electrostatic interaction between 
them, on the other hand, is simply 

Ug = i,j = +,~, (2) 

where £b is the Bjerrum length and takes the value of 
7 A in the systems henceforth considered. Here we will 
simply approximate the electrostatic interaction between 
the smeared-out charge distribution contained in the col- 
loidal disks as that between their centers of mass and 
thus assume a Coulomb form of interaction. Further, the 
configuration of the system is such that we can assume 
the three-dimensional form for this interaction. Refine- 
ment to the above assumption can be straightforwardly 
performed for other systems, such as those with salt. 

Since we only consider the mean-field limit in our 
model, the effects of electrostatics, which only comes in 
through the one-loop order, would not contribute to the 
gas phase. Therefore, the gas-phase free energy, / G , has 
two contributions: water-induced short-range repulsion 
and entropy. It takes the form 

2 

Pf G \k=o = Pilog^+p+log^+ +?-(4TrR 2 )w__ 
e e 2 

+ (4vra 2 ) w ++ + p + p_ {AiraR) t0-+, (3) 

where the colloidal number density p_ and counte- 
rion density p+ satisfy the electro-neutrality condition, 
Z + p + + Z-p- = 0. Again, the interaction between the 
water molecules and the colloidal particles can be en- 
coded directly into the parameters w,j characterizing the 
strength of these water-induced short-range repulsive in- 
teractions between the colloidal disks and counterions. 
Taking into account our conjecture, we simply assume 

that u>++ = u>__ < w |_. 

To demonstrate long-range attraction between nega- 
tively charged colloidal disks, we integrate out the fast 
modes associated with the smaller-sized counterions in 
order to obtain the effective interaction potential between 
the colloidal disks (l6j 

U—(k) = — + AttR 2 e' k R 

k 

- A Z_ (4) 

— t h 4ttci z w ++ e K a + -ri- 
fe P+ 

The structure factor S(k), in turn, is related to the 
effective interaction potential in Eq. 2] via S~ 1 (k) = 
U—{k) + pZ l . For a range of parameters, we do indeed 
observe the corresponding peak in the structure factor 



for this system. An example of which is shown in Fig. [2] 
Thus, we do see correlation over certain length scale be- 
tween the negatively charged colloidal disks in the gas 
phase of our system, a prelude to actual crystallization. 

For the periodic phase, it has been established that the 
optimal structure in the case of asymmetric charge ratios 
between the colloidal particles and counterions in two di- 
mensions is the hexagonal closed-packed lattice [10j. We 
model the two-dimensional closed-packed hexagonal lat- 
tice using a foam model developed earlier to study non- 
closed packed lattices in colloidal crystals [TJ E3- ^ n 
the foam model, the free energy for the solid is divided 
into an entropic contribution and an interaction contri- 
bution. To compute the entropic contribution, wc follow 
the method developed earlier [l3|, [3] by adopting the 
cellular free-volume theory. In this approximation, each 
colloidal particle is confined to a cage formed by its neigh- 
bors. The free volume available to each particle's center 
of mass is the volume of the Wigner-Seitz cell after a 
layer of thickness of a (where a is the hard-core radius of 
particles) is peeled off its faces. Using simple geometric 
arguments, the cellular free- volume approximation yields 
the following entropy contribution to the free energy in 
the high-density limit: 

Despite its mean-field nature, the free volume theory 
yields excellent quantitative agreement with available nu- 
merical simulations in the high-density limit [HI, [3] • 

On the other hand, the interaction contribution to the 
solid free energy can be carried out by direct summation 
of the reciprocal-lattice series [l5[ . To compute the inter- 
action energy per unit volume between the colloidal disks 
in the hexagonal lattice, we use the effective colloidal po- 
tential in Eq.[¥]and apply the Poisson summation formula 
that replaces the sum over the lattice vectors A by a sum 
over the reciprocal lattice vectors Q defined by the prop- 
erty Q ■ A = 2nm for some integer to. We thus obtain 

(3Uat = -^^*(Q) 2 C/-(-Q) , (6) 
lJK Q 

where A is the area of a unit cell and <r(fe) represents the 
transform of a circular-domain charge density 

m = 2 -^-Ji{kR) , (?) 

and Ji(.) is the regular Bessel function of order 1. In to- 
tal, the solid free energy is simply /3f s \ k=0 = Pf e + (3f int . 
Together with the gas free energy (3f \k=o in Eq. [31 we 
can determine phase coexistence of the periodic solid and 
dilute gas by using the common-tangent rule by imposing 
equal chemical potential and pressure for the two phases. 
The resulting phase diagrams are shown in Fig. [5] for par- 
ticle size-ratios of R/a = {6, 10, 14} and charge-ratios of 
\Z-/Z+\ = {10,15}. For all probed size-ratios, starting 
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FIG. 2: Phase Diagrams and Structure Factor. We show here the phase diagrams for size ratios of R/a = 14, R/a — 10, and 
R/a = 6 and for charge ratios of \Z\/Z2\ = 15 (blue dots) and \Z\/Zi\ = 10 (red circles). Starting from low colloidal number 
density p, we progress from the gas phase (G), to the coexistence region (G+S), to the solid phase (S). The close-packing 
number-density is p c — 0.2886/i? 2 (grey vertical line), which corresponds to the two-dimensional hexagonal close-packing 
volume fraction of <p = 0.9069. Here we assume that = W— = w. The structure factor is related to the effective 

interaction potential via S^ 1 (k) — U (fc) + pZ 1 in the low-density limit. This is done for a sample of colloidal disks with a 

size ratio of R/a = 14, a charge ratio of \Z—/Z+\ = 10, W-+/w = 6.8 and at a number density of p = 0.127/i? 2 within the 
coexistence region. 



at low colloidal number density p_ we progress from the 
gas phase (G), to the coexistence region (G+S), to finally 
the solid phase (S) at high density. The close-packing 
number-density is p c = 0.2886/i? 2 , which corresponds 
to the two-dimensional hexagonal close-packed volume 
fraction of f — 0.9069 (/ = nR 2 p). Here we assume that 

w++ = w = w, since our conjecture makes no a priori 

distinction between the two cases regarding their effects 
in water. 

From Fig. [5J we see that the coexistence region de- 
creases as the charge ratio increases from \Z-/Z+\ = 10 
to \Z-/Z+\ = 15. This is expected since when we in- 
crease the strength of the electrostatic attraction (Eq. 
[5]) between the negatively charged colloidal disks and 
the counterions, the water-induced attraction between 
the colloidal disks correspondingly decreases in relative 
strength, which allows for the increased gas region in the 
phase diagram. We also see that the same effect when we 
keep the electrostatic interaction constant and decrease 
the size-ratios between the colloidal disks and counteri- 
ons. To understand this, we recall that the short-range 
repulsion we model in Eq. ([T]) depends explicitly on the 



sizes of the particles. Therefore, decreasing their size- 
ratio is equivalent to decreasing the relative strength of 
their water-induced short-range repulsion. Hence, the 
same reason that swells the gas phase region in the phase 
diagram with increasing charge-ratio also applies here as 
we decrease the size-ratio. In this respect, our model 
provides concrete predictions that can be experimentally 
tested in its regime of applicability. 

Given the mean-field assumption in the construction 
of our gas and solid free-energy densities, our model does 
not indicate the existence of a critical point for the gas- 
solid transition in this compressible binary system involv- 
ing colloidal disks and counterions. At low enough ratios 

of w \./w (as a function of size- and charge-ratios), we 

see that the gas phase would abruptly preempt the solid 
phase by having a globally minimum free energy for all 
densities. The reverse is true in the opposite limit of high 

enough ratios of w 1- /w where the solid phase would 

preempt the gas phase. However, given that cellular 
free-volume approximation only works well in the high- 
density regime, we expect our model would no longer be 
valid in the phase-space region wherein the solid free en- 
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ergy starts dominating the gas free energy as the global 
minimum for all densities. For example, in the case of 
size-ratio of R/a — 14, the gas free energy becomes the 
global minimum for w \./w < 5, and our model indi- 
cates that the system would be completely gaseous. For 
the smaller charge-ratio of \Z-/Z+\ = 5, we find that 

our model becomes unreliable when w |_ /w > 10 for the 

size-ratio R/a = 6 at which point the solid free energy 
dominates the gas free energy for all densities, even in the 
low-density regime where the free- volume approximation 
is less reliable. Likewise, we approach the limit of appli- 
cability for our model when w y/w > 13 for R/a = 10, 

and w \-/w > 15 for R/a = 14. 

As discussed before, our model makes no a priori pre- 
dictions on the magnitudes of the water-induced short- 
range repulsion between different charge species. We can 
only infer the ordering of their magnitude. One possi- 
ble refinement of our mean-field model is to generalize 
the short-range repulsion parameters, w±± and w± T , as 
functions of particle size and molecular geometry. We 
can also account for the self-interaction between water 
molecules. Recently, there has been wide recognition in 
the important role water plays in many pertinent bio- 
logical processes [ijj ■ As a result, there has been much 
concerted effort on devising more accurate models for the 
various observed properties of water [l8j . 

We note that short-range attractions between 
macroions in the presence of high- valence counterions 
also lead to phase segregation [ljj]. These precipita- 
tion has been explained by various modified Poisson- 



Boltzmann approaches long ago [13, HH, HH, [H, HH . Our 
present model of segregation is very different in nature 
since it involves long-range attractions among negatively 
charged colloids due to assymetric hydration forces. Fur- 
ther, we have neglected charge renormalization here since 
ion condensation along the colloids has been found in di- 
lute three-dimensional colloidal suspensions [5^, [2(|; in 
our system, the lattice phase is concentrated. Moreover, 
charge renormalization would yield only negligible cor- 
rections to the theory. 

It is our hope that our simple model provides a con- 
crete example in which the molecular structure and 
chemical details of water contribute crucially to the un- 
derstanding of the phase properties of such systems as 
the self-assembly of colloidal disks. Having demonstrated 
another example of the inadequacy in treating water as 
a continuous, inert dielectric medium, it is our further 
hope that our model, with additional refinements, would 
serve as a starting point in modeling other more compli- 
cated, self-assembling systems such as charged colloidal 
suspensions in three dimensions [271 ] . 
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